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QJ ' Abstract 

The formation of singularities on a free surface of a conducting ideal fluid in a strong 
electric field is considered. It is found that the nonlinear equations of two-dimensional 
fluid motion can be solved in the small-angle approximation. This enables us to show 

P3 ■ that for almost arbitrary initial conditions the surface curvature becomes infinite in a 

t^l finite time. 

i 

C I Electrohydrodynamic instability of a free surface of a conducting fluid in an external 

Yl ' electric field [1,2] plays an essential role in a general problem of the electric strength. The 

•^ ' interaction of strong electric field with induced charges at the surface of the fluid (liquid 

>->. metal for applications) leads to the avalanche- like growth of surface perturbations and, as a 

'"5 ! consequence, to the formation of regions with high energy concentration which destruction 

can be accompanied by intensive emissive processes. 

In this Letter we will show that the nonlinear equations of motion of a conducting fluid 

^ ■ can be effectively solved in the approximation of small perturbations of the boundary. This 

allows us to study the nonlinear dynamics of the electrohydrodynamic instability and, in 



O 

^^ . 

O . particular, the most physically meaningful singular solutions. 

O 



Let us consider an irrotational motion of a conducting ideal fluid with a free surface, 
O I ^ = vi^^y^t)^ that occupies the region — oo < z < ri{x,y,t), in an external uniform electric 

^ ■ fleld E. We will assume the influence of gravitational and capillary forces to be negligibly 

5j . small, which corresponds to the condition 
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where g is the acceleration of gravity, a is the surface tension coefficient, and p is the mass 
density. 
^ ■ The potential of the electric fleld cp satisfles the Laplace equation. 



with the following boundary conditions, 

(f — > —Ez, 

f = 0, Z = 7]. 



The velocity potential $ satisfies the incompressibility equation 

A$ = 0, 
which one should solve together with the dynamic and kinematic relations on the free surface, 

9$ (V$)2 (Vv?)2 
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z = V, 



dr] (9$ 
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z = ri, 



dt dz 
where F is some function of variable t, and the boundary condition 



$-^0, 



-oo. 



The quantities 77(0;, y, t) ipi^x, y, t) = ^\z=ri are canonically conjugated, so that the equations 
of motion take the Hamiltonian form [3], 
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where the Hamiltonian 
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coincides with the total energy of a system. With the help of the Green formula it can be 
rewritten as the surface integral. 
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where (f = ip + Ez is the perturbation of the electric field potential; ds is the surface 
differential. 

Let us assume {Vrjl <^ 1, which corresponds to the approximation of small surface angles. 
In such a case we can expand the integrand in a power series of canonical variables rj and ip. 
Restricting ourselves to quadratic and cubic terms we find after scale transformations 

t^tE-\ATrpy/^, ^^^E/{ATrpy/^, H ^ HE^/{ATrp) 

the following expression for the Hamiltonian, 
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Here k is the integral operator with the difference kernel, whose Fourier transform is the 
modulus of the wave vector, 

+00+00 , 
kf = -—f I i^^^^^ :-^dx'dy'. 

The equations of motion, corresponding to this Hamiltonian, take the following form, 

A-kri = - [{ki'f - (V^)^ + {krif - (Vrif] + k{rikri) + V{ri\/ri), (1) 

r]t-kilj = -k{'r]kilj)-V{'r]Vip). (2) 

Subtraction of Eqs. (2) and (1) gives in the linear approximation the relaxation equation 

{ip-r])t = -kC^-v), 

whence it follows that we can set ip = rj in the nonlinear terms of Eqs. (1) and (2), which 
allows us to simplify the equations of motion. Actually, adding Eqs. (1) and (2) we obtain 
an equation for a new function f = {ip + ri)/2, 

ft-kf = \ {kff - \ (V/)^ (3) 

which corresponds to the consideration of the growing branch of the solutions. As / = r^ in 
the linear approximation, Eq. (3) governs the behavior of the elevation 77. 

First we consider the one-dimensional case when function / depends only on x (and t) 
and the integral operator k can be expressed in terms of the Hilbert transform H ^ 
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where P denotes the principal value of the integral. As a result, Eq. (3) can be rewritten as 

ft + HU = \{HUf-\{Uf. (4) 

It should be noted that if one introduces a new function / = Hf, then Eq. (4) transforms 
into the equation proposed in Ref. [4] for the description of the nonlinear stages of the 
Kelvin-Helmholtz instability. 

For further consideration it is convenient to introduce a function, analytically extendable 
into the upper half-plane of the complex variable x. 

Then Eq. (4) takes the form 

Re {vt + ivx + 2f f J.) = 0, 



that is, the investigation of integro-differential equation (4) amounts to the analysis of the 
partial differential equation 

vt + iv^ + 2VV.J, = 0, (5) 

which describes the wave breaking in the complex plane. Let us study this process in analogy 
with [5,6], where a similar problem was considered. Eq. (5) can be solved by the standard 
method of characteristics, 

v = Qix'), (6) 

x = x' + it + 2Q{x')t. (7) 

where the function Q is defined from initial conditions. It is clear that in order to obtain an 
explicit form of the solution we must resolve Eq. (7) with respect to x'. A mapping x -^ x', 
defined by Eq. (7), will be ambiguous if dx/dx' = in some point, i.e. 

1 + 2Q^d = 0. (8) 

Solution of (8) gives a trajectory x' = x'{t) on the complex plane x'. Then the motion of the 
branch points of the function v is defined by an expression 

x{t) = x'{t) +it + 2Q{x'{t))t. 

At some moment to when the branch point touches the real axis, the analiticity of v{x,t) 
at the upper half-plane of variable x breaks, and a singularity appears in the solution of 
Eq. (4). 

Let us consider the solution behavior close to the singularity. Expansion of (6) and (7) 
at a small vicinity of x = x(to) up to the leading orders gives 

v = Qo- 6x'/i2to), 

Sx = i6t + 2Qo6t + Q"to{6x')^, 

where Qo = Qi^'ito)), Q" = Qx'x'{x'(to)), 6x = x — x(tQ), 6x' = x' — x' (to) , and 6t = t — tQ. 
Eliminating 6x' from these equations, we find that close to singularity v^ can be represented 
in the self-similar form {6x ~ 6t), 
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As Re(f ) = ri/2 in the linear approximation, we have at t = to 

Vxx ~ \Sx\'^^^, 

that is the surface curvature becomes infinite in a finite time. It should be mentioned that 
such a behavior of the charged surface is similar to the behavior of a free surface of an ideal 
fluid in the absence of external forces [5,6], though the singularities are of a different nature 
(in the latter case the singularity formation is connected with inertial forces). 



Let us show that the solutions corresponding to the root singularity regime are consistent 
with the applicability condition of the truncated equation (3). Let Q{x') be a rational 
function with one pole in the lower half-plane, 

Qi^') = - ,, ''.,,2 > (9) 

2[x + lAy 

which corresponds to the spatially localized one-dimensional perturbation of the surface 
(s > and A > 0). The characteristic surface angles are thought to be small, 7 ~ s/A^ <^ 1. 
It is clear from the symmetries of (9) that the most rapid branch point touches the real 
axis at X = 0. Then the critical moment to can be found directly from Eqs. (7) and (8). 
Expansion of to with respect to the small parameter 7 gives 

to ^ A [1 - 3(7/4)1/=^] . (10) 

Taking into account that the evolution of the surface perturbation can be described by an 
approximate formula 

s{A-t) 

we have for the dynamics of the characteristic angles 
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Then, substituting the expression for to (10) into this formula, we find that at the moment 
of the singularity formation with the required accuracy 

7(to) ~ 7'/', 

that is, the angles remain small and the root singularities are consistent with our assumption 
about small surface angles. 

In conclusion, we would like to consider the more general case where the weak dependence 
of all quantities from the spatial variable y is taken into account. One can find that if the 
condition \kx\ <^ \ky\ holds for the characteristic wave numbers, then the evolution of the 
fiuid surface is described by an equation 

[Vt + iVx + 2vVa:]^ = -iVyy/2, 

which extends Eq. (5) to the two-dimensional case. 

An interesting group of particular solutions of this equation can be found with the help 
of substitution v{x, y, t) = w{z, t), where 

i (y-yo? 

Z = X — - 



2 t 
The equation for w looks like 

Wt + iwz + 2wWz = —w/{2t). 
5 



It is integrable by the method of characteristics, so that we can study the analyticity violation 
similarly to the one-dimensional case. Considering a motion of branch points in the complex 
plane of the variable z we find that a singularity arises at some moment to < at the point 
I/O along the y-axis. Close to the singular point at the critical moment t = to we get 

1/2 ^1 _ |A„|-l 



Vxx\sy=o ~ \0X\ , Vxx\sx=0 

This means that in the examined quasi-two-dimensional case the second derivative of the 
surface profile becomes infinite at a single isolated point. 

Thus, the consideration of the behavior of a conducting fluid surface in a strong electric 
field shows that the nonlinearity determines the tendency for the formation of singularities 
of the root character, corresponding to the surface points with infinite curvature. We can 
assume that such weak singularities serve as the origin of the more powerful singularities 
observed in the experiments [7,8]. 
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